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1. Introduction 

Transport of particles and energy across the confining magnetic field of fusion devices 
is anomalous [Q, i.e., it is much larger than the neoclassical transport level associated 
with binary collisions in a toroidal geometry [2j. It is thought that anomalous transport 
is caused by plasma turbulence, which in turn manifests itself as fluctuations in most 
plasma parameters. To understand anomalous transport, a two-pronged approach is 
being applied: (i) sophisticated diagnostics measure fluctuations and (ii) advanced 
simulations are being developed and compared to these measurements. Once our 
understanding of the relationship between fluctuation-induced anomalous transport and 
plasma confinement quality is more complete, we will be able to reduce transport due 
to the identified mechanism (s). 

The fusion plasma measurements presented in this paper are of fluctuations in the 
electron density. Small-angle collective scattering [3J Hj was used in the Wendelstein 
7-AS (W7-AS) stellarator [5 j and phase-contrast imaging (PCI) |E] is being used in the 
Alcator C-Mod tokamak [7j. 

We specifically study density fluctuation power versus wavenumber (also known 
as the wavenumber spectrum) in W7-AS and C-Mod. These wavenumber spectra 
characterize the nonlinear interaction between turbulent modes having different length 
scales. 

The second part of our measurements, wavenumber spectra (i) of galaxies from the 
Sloan Digital Sky Survey (SDSS) 8j and (ii) from a variety of sources (including the 
SDSS data) are published in Ref. P and have been made available to us [TP] . 

The paper is organized as follows: In section El we review our initial results from 
Ref. fl] . Thereafter we analyze our expanded data set in section El and in response to 
the results revise our treatment of the original W7-AS measurements in section 0] A 
discussion follows in section El and we conclude in section |H1 

2. Review of earlier findings 

2.1. Stellarator fusion plasmas 

In Ref. jTl] we proposed that the density fluctuation power P decreases exponentially 
with increasing wavenumber k on mm scales in fusion plasmas according to 



where n > is a constant having a dimension of length and k = 2tt/\, where A is the 
corresponding wavelength. Initially we suggested the simplified form 



P(k) oc - x e 



—nk 



(1) 



P(k) oc e 



—nk 



(2) 



in Ref. [12]. Eq. (J2J) also holds for density fluctuations in the Tore Supra tokamak |13j . 
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A wavenumber spectrum of turbulence in W7-AS is shown in Fig. ^ The measured 
points are shown along with two power-law fits 

P{k) oc k~ m , (3) 

where m is a dimensionless constant. The power-law fits are shown as solid lines and 
an exponential fit to Eq. is shown as a dashed line. 

All fits shown in this paper have a normalized x 2 < 1, ensuring the quality of the 
fits is good. The error bars are standard deviations and the semi-transparent rectangles 
indicate which points are included to make the fits. 

The power-law fits are motivated by classical fluid turbulence theory where one 
expects wavenumber spectra to exhibit power-law behavior with exponents m depending 
on the dimension of the observed turbulence: 

• 3D: Energy is injected at a large scale and redistributed (cascaded) by nonlinear 
interactions down to a small dissipative scale. In this case, the energy spectrum in 
the inertial range E{k) oc k~ 5 ^ 3 |14j . 

• 2D: Here, two power-laws exist on either side of the energy injection scale. For 
smaller wavenumbers, the inverse energy cascade obeys E{k) oc k~ 5 ^ 3 and for larger 
wavenumbers, the enstrophy cascade follows E(k) oc k~ 3 [To] . 

• ID: Energy is injected at a large scale and dissipated at a small scale; E{k) oc k~ 2 

m- 

Our measured power is equivalent to the d- dimensional energy spectrum F d {k) 

pud US] 

P(k) = F d [k) = ^ 

A 1 = 2 A 2 = 2nk A 3 = 4vrfc 2 , (4) 

where A d is the surface area of a sphere having radius k and dimension d. Usually on 
would assume that d = 2 in fusion plasmas, since transport along magnetic field lines is 
nearly instantaneous. The fits to Eq. (jH)) in Fig. ^yield exponents m = 3.0 ± 0.4 (small 
wavenumbers) and 6.9 ± 0.7 (large wavenumbers). A similar behavior has previously 
been reported in Ref. jTH] where it was speculated that the wavenumber value at the 
transition between the two power-laws should correspond to a characteristic spatial scale 
in the plasma. The only length scale close to the transitional value was found to be the 
ion Larmor radius pi. 

The spectrum at small wavenumbers is roughly consistent with the inverse energy 
cascade in 2D turbulence, i 7 ^^) oc k~ 8 ^ 3 . The exponent at large wavenumbers does not 
fit into this framework. However, for very large wavenumbers one enters the dissipation 
range; here, it has been argued that the energy spectrum could have one of the following 
dependencies 

-ENeumann(^) OC C -^Heisenberg (k) OC k 7 , (5) 



Density fluctuations on mm and Mpc scales 



4 



where a > is a constant having a dimension of length (see Ref. ^B] and references 
therein). The energy spectrum proposed by J von Neumann was what initially inspired 
us to investigate an exponential decay of P(k) in Ref. ^2j- J von Neumann's work on 
this topic is from 1949 and two years later, A A Townsend published a more generalized 
expression j2Hj. 

Fitting all wavenumbers to Eq. (JTJ), E Neumarin (k) / A 2 , we find that n = 0.11 ± 0.004 
cm or a wavenumber of 57.1 cm -1 . Alternatively, the transitional wavenumber found 
at the intersection between the two power-laws is 32.8 cm -1 (0.19 cm). The expression 
En eisenherg (k) / A 2 yields m — 8, which is close to the experimental value m = 6.9 for 
large wavenumbers. Calculating the ion Larmor radius at the electron temperature p s 
for this case we find that it is 0.1 cm, i.e. the same order of magnitude as the spatial 
scales found above. We used p s instead of pi because ion temperature measurements 
were unavailable. 

Currently we can think of three possible explanations for the behavior of the 
wavenumber spectrum: 

(i) We observe 2D turbulence and the transition between the two power-laws occurs at 
a spatial scale where the inverse energy cascade develops into the dissipation range. 
However, the enstrophy cascade is not accounted for in this case. 

(ii) We observe 2D turbulence in the dissipation range described by Eq. (JTJ) as proposed 
by J von Neumann. 

(iii) Turbulence theory does not apply. The transition between two power-laws or the 
characteristic scale found using a single exponential function (Eq. ©) indicates 
that one scale dominates the turbulent dynamics in the wavenumber range studied. 

In measurements of 2D fluid turbulence, it has been demonstrated that the inverse 
energy and forward enstrophy cascades merge into a single power-law when the system 
transitions to being fully turbulent [21]. This might be the reason for the missing 
enstrophy cascade discussed in item (i) above. 

2.2. Cosmology 

The shape of the wavenumber spectrum shown in Fig. ^ bears a striking resemblance 
to measurements of fluctuations in the distribution of galaxies presented in Ref. , see 
Fig. |21 This motivates us to apply the analysis described in section 12.11 to the galaxy 
data. In section 12.2.11 we briefly put these measurements in context and then present 
an analysis of the galaxy wavenumber spectrum in section 12.2.21 

2.2.1. Inflation Dramatic developments have taken place in cosmology over the last 
decade, lending increasing support to the paradigm of inflation as an explanation for 
what took place before the events described by the big bang theory [22] ■ Inflation solved 
the so-called horizon and flatness problems, but was at odds with earlier observations 
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indicating that the ratio of the mass density of the universe to the critical value, the 
density parameter Q, was 0.2-0.3, while inflation predicted it should be 1: 

Pc 

Q < 1 : open Q — 1 : flat Q > 1 : closed, (6) 

where p c = 3Hq/8ttG is the critical mass density, H ~ 70 km/s/Mpc is the Hubble 
parameter observed today and G is I Newton's gravitational constant • However, new 
measurements in the late 1990's lead to a drastic modification of Q: Observations of type 
la supernovae (SN) showed that the separation velocity between galaxies was speeding 
up, not slowing down as would be expected for an open universe. The underlying 
explanation for this accelerated expansion is not understood, but it seems that the 
universe contains large quantities of negative pressure substance, creating a gravitational 
repulsion driving the expansion. This negative pressure material is called dark energy, 
the total density of dark energy Q\ is 0.7. The existence of dark energy is equivalent 
to the cosmological constant A introduced by A Einstein. The dark matter density 
Qd is 0.25 and the baryonic matter density is 0.05, so the total density is very 
close (or equal) to the critical density. The SN la data is supported by measurements 
of nonuniformities in the cosmic microwave background (CMB) radiation. The CMB 
anisotropy is due to the presence of tiny primordial density fluctuations at the time of 
recombination, where atoms formed. At that point in time the age of the universe was 
about 400,000 years and the temperature was 3000 K. The structures observed in the 
CMB are called acoustic peaks, and the simplest versions of inflation all reproduce these 
structures quite accurately. The acoustic peaks can not be modelled by assuming that 
the universe is open. 

2.2.2. SDSS wavenumber spectrum A study of density fluctuations on large scales 
using 205,443 galaxies has been published by the SDSS Team in Ref. f"||, see Fig. [21 
3D maps of the universe are provided by the SDSS galaxy redshift survey, observing 
about a quarter of the celestial sphere using a 2.5 m telescope and a charge-coupled 
device (CCD) camera. The galaxies had a mean redshift z ~ 0.1, corresponding to 
light emitted 1-2 Gyr ago Fixing some cosmological parameters to Wilkinson 

Microwave Anisotropy Probe (WMAP) satellite values [21] one finds - using physics 
based models - that the wavenumber spectrum measurements were fitted by a matter 
density Q m = Q d +Q b = 0.295±0.0323. In this case h = # /(100 km/s/Mpc) = 0.72 was 
assumed and it was observed that the wavenumber spectrum was not well characterized 
by a single power-law. 

A follow-up paper by the SDSS Team, Ref. [2*5] . combined non-CMB measurements 
(SDSS) with CMB measurements (WMAP) to constrain free parameters of cosmological 
models and break CMB degeneracies in parameter space. This resulted in Q m = 
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0.30 ±0.04 and h = 0.70±g;g|. Adding the SDSS information more than halved WMAP- 
only error bars on some parameters, e.g. the Hubble parameter and matter density. 

The data presented in Fig. [21 has been taken from M Tegmark's homepage |2fij . 
According to the recommendation by the SDSS Team jjjj, the three largest wavenumbers 
shown are not used in the fits described in this section. 

As we did for the W7-AS data in Section 12.11 we fit the SDSS measurements to 
two power-laws (Eq. Q) or the exponential function with an algebraic pre-factor (Eq. 
(0)). The power- law fits are shown as solid lines, the exponential fit as a dashed line. 

The power-law fits yield exponents m = 0.8 ± 0.03 (small wavenumbers) and 1.4 ± 
0.1 (large wavenumbers). The wavenumber ranges were determined by minimizing the 
combined x 2 of the fits. As the SDSS Team found, a single power-law can not describe 
the observations. The exponents are not close to the ones governing fluid turbulence 
discussed in section l2~T1 The transitional wavenumber (power-law intersection) is 0.07 
h Mpc -1 , corresponding to a length of 89.8 h -1 Mpc. 

We find the characteristic length from Eq. (0) to be n = 2.3 ± 0.5 h -1 Mpc or a 
wavenumber of 2.7 h Mpc -1 . 

3. Analysis of additional measurements 

In Ref. [II] w e noted that both fusion plasma and cosmological wavenumber spectra 
peak at small wavenumbers and decrease both above and below that peak. In this section 
we analyze supplemental measurements on a wider range of scales. The fusion plasma 
measurements were made in a tokamak, so our explicit assumption is that turbulence 
in stellarators and tokamaks is comparable. 

We fit the data to a modified version of Eq. flj) . namely 

P{k) cx -I x c -* (7) 

where a is brought in as an additional fit parameter. The introduction of a is based 
on the assumed functional form of the energy spectrum in the dissipation range of fluid 
turbulence [2*7] . Further, since we have no clear criteria for choosing between Eqs. (Q) 
and (j2J), leaving a to be fitted will test our biased opinion in Ref. [TT] for Eq. (JTJ). Using 
Eq. ®, Eq. © implies that E(k) oc k^ 1 ^ x e~ nk . 

3.1. Tokamak fusion plasma 

We analyze density fluctuations at small wavenumbers in C-Mod using the PCI 
diagnostic. Measurements in the low confinement mode taken from Fig. 11 in Ref. 
[2H1 are shown in Fig. 01 The dashed line shows the fit made using Eq. (JZj) which yields 
n = 0.11 ± 0.005 cm and a = 0.30 ± 0.02. The size of n is close to the value of p s in 
the given plasma. 
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3.2. Cosmological wavenumber spectrum 

The expanded cosmological data set shown in Fig. 38 of Ref. is re-plotted in Fig. |U 
These measurements have been provided by M Tegmark ^0] • They are a combination of 
the measured density fluctuation power using several diagnostic techniques; the dashed 
line is a fit to a current model in cosmology, the so-called "vanilla light" flat adiabatic 
ACDM model with negligible tilt, running tilt, tensor modes or massive neutrinos |9~ll25|. 

During our analysis of the measurements it became obvious that Eq. (jJJ) does not 
fit the entire range of wavenumbers above the peak of the spectrum. Therefore two fits 
were made: 

• A small wavenumber fit where the smallest wavenumber used is chosen to be the 
one closest to (but larger than) the peak of the ACDM model shown in Fig. |3J The 
largest wavenumber used is chosen so that \ 2 of the fit is 1 . In Fig. the fit is 
displayed as a dashed line. The parameters found are n = 13.3 ± 1.6 h -1 Mpc and 
a = 0.24 ± 0.06. Note the interesting proximity of this a to the result in section 

o 

• A large wavenumber fit where the largest wavenumber used is the maximum of the 
data set and the smallest wavenumber is chosen so that x 2 °f the fit is 1. In Fig. 
El the fit is displayed as a dashed line. The parameters found are n = 0.33 ± 0.05 
h- 1 Mpc and a = 1.90 ± 0.05. 

The wavenumbers not being used in either fit define a transitional region: [0.08, 
0.10] h Mpc -1 . This interval is in rough agreement with the transition found between 
the two power-law fits shown in Fig. 

4. Revised analysis of stellarator fusion plasmas 

Based on the insights gained in section El we revisit the turbulence measurements made 
in W7-AS, see Fig. ^ Our approach is to make three fits, one for small, one for large 
and one for all wavenumbers: 

• For the fit to small wavenumbers we use n and a found from the PCI measurements 
in section l3~T1 and only allow the overall amplitude to vary. The largest wavenumber 
to be included in the fit is chosen so that x 2 °f the fit is 1. Using this criterion the 
four smallest wavenumbers are utilized. In Fig. 0the fit is shown as a dashed line. 

• We allow all fit parameters in Eq. (J2J) to vary for the fit to the four largest 
wavenumbers and show the resulting fit as a dashed line in Fig. |H1 The fit yields 
n = 0.14 ± 0.008 cm and a = 0.25 ± 0.1. This fit is also good for the entire 
wavenumber range: Using the found parameters, \ 2 remains less than 1 for all 
eight wavenumbers. 

• Fit to Eq. (JJJ) using all eight wavenumbers and allowing every fit parameter to 
vary. The fit yields n = 0.14 ± 0.007 cm and a = 0.11 ± 0.2, which is consistent 
with a pure exponential function (Eq. ((2))). 
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The amplitude-only fit using PCI parameters indicates a transition at the center 
of the data set, whereas the two other fits suggest that the data is consistent with a 
continuous behavior. It is apparent that additional measurements at large wavenumbers 
in fusion plasmas are needed to gain confidence in the fits. Some of the existing points 
might be in a transitional region akin to the one found in the cosmological data. 

5. Discussion 

The fact that density fluctuations on small (fusion plasmas) and large (cosmological) 
scales can be described by an exponential function with an algebraic pre-factor, 
Eq. (J2J), might indicate that plasma turbulence at early times has been expanded 
to cosmological proportions. A natural consequence of that thought would be to 
investigate fluctuations in quark-gluon plasmas (QGPs) corresponding to even earlier 
times. However, experimental techniques to do this are not sufficiently developed at the 
moment due to the extreme nature of QGPs. It has been suggested that complex (or 
dusty) plasmas could be used as a model of QGPs [29J. 

It is fascinating that wavenumber spectra over wider scales peak at small 
wavenumbers and decrease both above and below that peak. This is seen both in 
fusion plasmas and on cosmological scales, compare Figs. HU and 0J Turbulence theory 
in ID or 3D would interpret the peak position as the scale where energy is injected. 

Fitting wavenumber spectra to power-laws as we did in section El is based on fluid 
turbulence theories, but in general care must be taken when interpreting the outcome: 
We know that an exponential function can be Taylor expanded to a Maclaurin series 



So locally, i.e. for a small range of wavenumbers, an exponential dependency can 
be masked as a power-law; the exponent would vary as a function of the wavenumber 
range selected. To some extent this objection also holds true for our results using 
Eq. ((7j) to identify two wavenumber intervals having different fit parameters. It is 
an open question whether an extension of the cosmological data to larger wavenumbers 
would necessitate further intervals. However, in fluid turbulence simulations a transition 
between near- and far-dissipation similar to the one we discuss has been identified [30J. 
This observation lends support to the interpretation that our dual fits describe different 
types of dissipation. 

To sum up, we favor Eq. as a descriptor of the data, both for fusion plasma and 
cosmological measurements. Perhaps forcing occurs at a large scale where the spectra 
peak and transitions directly to dissipation. The fact that we obtain n ~ p s using Eq. 
(J7J) for fusion plasmas suggests that n is the characteristic scale of turbulence. 

We would like to point out that observations of electric field fluctuations in the solar 
wind are found to be proportional to exp(— kpi/12.5) for kpi > 2.5 |HT| . It should be 
noted that the electric field measurements at these large wavenumbers are noisy. Work 
on the solar wind data to analyze density fluctuations is in progress |32j- 




(8) 
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6. Conclusions 

We have in this paper reported on suggestive similarities between density fluctuation 
power versus wavenumber on small (mm) and large (Mpc) scales. 

The small scale measurements were made in fusion plasmas and compared to 
predictions from turbulence theory. The data sets fit Eq. (jJJ), which has a functional 
form that can be explained as dissipation by turbulence theory. The large scale 
cosmological measurements can also be described by Eq. (JZj). In general, two 
wavenumber ranges separated by a transitional region are identified. 

The similar dependency of density fluctuation power on wavenumber might indicate 
a common origin of these fluctuations, perhaps from fluctuations in QGPs at early stages 
in the formation of the universe. The value of a is almost identical for both fusion 
plasma and cosmological measurements at wavenumbers close to but above the peak of 
the spectra. 

To progress further, it is essential that the quantity of wavenumber-resolved fusion 
plasma turbulence measurements is vastly increased. 
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Figure 1. Wavcnumber spectrum of broadband turbulence in W7-AS. Squares are 
measured points. Solid lines are power-law fits (Eq. and the dashed line is a fit to 
Eq. The power-law fit grouping is indicated by the semi-transparent rectangles, 
the fit to Eq. Q uses all data points. 
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Figure 2. Wavenumber spectrum of galaxies measured by the SDSS Team. Squares 
are measured points. Solid lines are power-law fits (Eq. {2Jl), the dashed line is a fit 
to Eq. The power-law fit grouping of points is chosen so the combined \ 2 °f 

the fits is minimized. The power-law fit grouping is indicated by the semi-transparent 
rectangles, the fit to Eq. uses all data points up to k = 0.2 h Mpc^ 1 . The data set 
is taken from Ref. 
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Figure 3. Wavenumber spectrum of broadband turbulence in C-Mod. Squares are 
measured points. The dashed line is a fit to Eq. Q using the measurements indicated 
by the semi-transparent rectangle. 
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Figure 4. Wavenumber spectrum of the combined cosmological measurements. 
Asterisks are CMB radiation, squares are SDSS, plus is cluster abundance, diamonds 
are weak lensing and triangles are Lyman-a forest measurements. The dashed line is 
a fit to the ACDM model. The data sets are taken from Ref. [Sj. 
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Figure 5. Wavenumber spectrum of the combined cosmological measurements. 
Squares are measured points. The dashed line is a fit to Eq. (J7J using the data 
indicated by the semi-transparent rectangle. The measurements are taken from Rcf. 

Q. 
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Figure 6. Wavenumber spectrum of the combined cosmological measurements. 
Squares are measured points. The dashed line is a fit to Eq. @ using the data 
indicated by the semi-transparent rectangle. The measurements are taken from Rcf. 

Q. 
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Figure 7. Wavenumber spectrum of broadband turbulence in W7-AS. Squares are 
measured points. The dashed line is a fit to Eq. Q using the measurements indicated 
by the semi-transparent rectangle. Note that this is exclusively an amplitude fit, we 
use n and a from the PCI fit, see section IXT1 and Fig. |3] 
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Figure 8. Wavcnumber spectrum of broadband turbulence in W7-AS. Squares are 
measured points. The dashed line is a fit to Eq. Q using the measurements indicated 
by the semi-transparent rectangle. 



